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ABSTRACT 



Context. The ESO-Spitzer extragalactic Imaging Survey (ESIS) is the optical follow up of the Spitzer Wide-Area InfraRed 
Extragalactic (SWIRE) survey in the ELAIS-Sl area. 

Aims. The multiwavelength study of galaxy emission is the key to understand the interplay of the various components of 
galaxies and to trace their role in cosmic evolution. ESIS provides optical identification and colors of Spitzer IR galaxies and 
builds the bases for photometric redshift estimates. 

Methods. This paper presents B, V, R Wide Field Imager observations of the first 1.5 square degree of the ESIS survey. 
Data reduction is described including astrometric calibration, illumination and color corrections. Synthetic sources are 
simulated in scientific and super-sky-fiat images, with the purpose of estimating completeness and photometric accuracy for 
the survey. Number counts and color distributions are compared to literature observational and theoretical data, including 
non-evolutionary, PLE, evolutionary and semi-analytic ACDM galaxy models, as well as Milky Way stellar predictions. The 
ELAIS-Sl area benefits from extensive follow-up from X-ray to radio frequencies: some potential uses of the multi-wavelength 
observations are illustrated. 

Results. Object coordinates are defined with an accuracy as good as ~ 0.15 [arcsec] r.m.s. with respect to GSC 2.2; flux uncer- 
tainties are ~2, 10, 20% at mag. 20, 23, 24 respectively (Vega); we reach 95% completeness at B,V~25 and R~24.5. ESIS galaxy 
number counts are in good agreement with previous works and are best reproduced by evolutionary and hierarchical ACDM 
scenarios. Optical-Spitzer color-color plots promise to be very powerful tools to disentangle different classes of sources (e.g. 
AGNs, starbursts, quiescent galaxies). Ultraviolet GALEX data are matched to optical and Spitzer samples, leading to a discus- 
sion of galaxy properties in the UV-to-24 /im color space. The spectral energy distribution of a few objects, from the X-rays to 
the far-IR are presented as examples of the multi- wavelength study of galaxy emission components in different spectral domains. 



Key words. Surveys - Galaxies: evolution - Cosmology: observations - Infrared: galaxies - Ultraviolet: galaxies - Galaxies: 
statistics 
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* Based on observations collected at the European Southern 
Observatory, Chile, ESO No. 168 .A-0322(A). 

** ESIS web page: |http://dipastro.pd.astro.it/esis/| 



1. Introduction 

The assembly and evolution of galaxies, galaxy clusters 
and cosmic large-scale structure (LSS) are currently ma- 
jor issues in both theoretical and observational cosmol- 
ogy. Tracing the history of cosmic star formation and the 
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growth of the cosmic stellar mass density will lead to an 
understanding of the fundamental processes transforming 
the primordial diffuse plasma into the highly structured 
present-day Universe. 

The multi- wavelength study of the emission of galaxies 
provides tools to disentangle their different physical com- 
ponents. Young stellar populations power the UV-optical 
light, old stars emit predominantly in the near-IR, while 
dust — heated by either starburst activity or an active 
galactic nucleus (AGN) — dominates the mid- and far-IR 
luminosity. X-ray and radio luminosities are also produced 
by starbursts or AGNs. Extending the analysis of galaxy 
properties to as wide a redshift range as possible is com- 
pulsory, in order to understand the interplay of the various 
components and to trace their role in cosmic galaxy evo- 
lution. 

The two main instruments onboard Spitzer, observing 
in the mid- to far-IR (3 - 8 /^m for IRAC and 24, 70, 160 
jum for MIPS), were specifically designed to probe the old 
stellar content and dust re-radiation from distant galaxies. 

T he cosmic Infra- Red Background fCIRB. IPuget et aP 

llflflfitlHaniser et a,1.lll discovered by the COBE satel- 
lite in the mid 90's. is the most en erget ic diffuse radiation 
after t he CMB. lElbaz et alJ l)2002l) and lPranceschini et alJ 
lj2003|) have shown that at least 50% of the CIRB emis- 
sion is powered by luminous and ultra-luminous massive 
star-forming galaxies, at redshift z = 0.5 — 1.5, strongly 
evolving in cosmic time. The Universe seems to have ex- 
perienced a phase of enhanced activity of star-formation 
and gravitational accretion in the past, m ostly visible in 
the infrared fe.g. iFranceschini et al.lf200l[l . 

The Spitzer M ultiband Imaging Photometer (MIPS, 
iRieke et al.l 120041) provides the observations needed to 
push the study of mid- and far-IR sources to fainter 
luminosities and larger distances. Polycyclic Aromatic 
Hydrocarbon (PAH) emission, typical of starburst 7-13 
/zm restframe spectra, is sampled by MIPS up to redshift 

3, corresponding to an epoch when the Universe had 
only ~15% of today's age. 

A central issue in modern cosmology is when galax- 
ies assemble d their baryonic mass. The h ierarchical sce- 
nario (e.g. iKaiiffmann fc CharlotI Il998l) predicts that 
the most massive systems (e.g. Mgi^^j-g > 10^^ M©) 
formed relatively late through a slow process of merg- 
ing of smaller ga laxies. In monolithic-collapse models 
l|Eggen et alJll962l) . the bulk of stars formed in early-type 
galaxies at very high redshift, while subsequent merging 
and star formation are limited. Mo dern ACDM mo dels 
and hydrodynamical simulations (e.g. lNagaminel200l|) im- 
plement a mixture of the two. 

Several studies have attempted to investigate the for- 
mation and evolution of massive systems an d test model 
predictions (see, e.g.. the recent works bv JProrv et alJ 
l2005llTreul2nn4lFonta,na et a,].l2004ICima,ttil200,'^. but a 
clear and coherent picture has not yet emerged. Very little 
is known ab out galaxy stellar mass assembly at redshift 
z > 1.5 (e.g. lCimattill2003HDaddi et al.ll2004l) . 



The Infrared Array Camera (IRAG. IPazio et al.ll20o4) 
onboard Spitzer observes in the 3.6 — 8.0 fim wavelength 
range. The instrument was specifically designed for de- 
tecting galaxies' restframe near-IR emission, up to redshift 
z > 3, hence directly probing their stellar mass assembly. 

1.1. The SWIRE & ESIS Surveys 

A significant fraction of the first year of Spitzer in-flight 
operations has been devoted to six different Legacy science 
Programs, representing projects of general and lasting im- 
portance to the broad astronomical community. Among 
these. SWIRE (.Lonsdale et al...200,'l. .20041 and GOODS 
l|Dickinson et al.ll2003() are dedicated to cosmology. 

The Spitzer Wide-area Infra-Red Extragalactic survey 
(SWIRE^) is the largest Spitzer Legacy Program. It con- 
sists of a wide-area, imaging campaign designed to trace 
the evolution of dusty, star-forming galaxies, evolved stel- 
lar populations, and AGN as a function of environment, 
from redshifts z ~ 3, down to the current epoch. SWIRE 
includes 6 high-latitude fields, totaling 49 [deg^] in all the 
seven Spitzer bands. 

The large sky area covered by the SWIRE survey im- 
plies large galaxy samples of all kinds, secures high statis- 
tical significance to clustering and galaxy evolution stud- 
ies, and allows the possibility of rare object searches. 

The ESO-Spitzer wide-area Imaging Survey (ESIS) is 
an ESO Large Programme (P.I. Alberto Franceschini), 
securing optical ground-based imaging follow-up to the 
SWIRE Spitzer survey in the ELAIS-Sl field. 

The ELAIS-Sl region, together with Lockman Hole, is 
the highest priority region for SWI RE, thanks to the ex - 
tremely low 100 /im cirrus emission (ISchlegel et al ."l998V 
In fact, it includes the absolute minimum of the Galactic 
100 emission in the Southern sky. 

A broad wavelength coverage is required by the nature 
of the sources we are targeting, which are typically red- 
dened, and/or with old evolved stellar populations, and 
at high redshifts. These properties imply dramatic diffi- 
culties when trying to acquire spectroscopic follow-up in 
the optical. Photometric investigation of the properties of 
faint sources provides a very powerful alternative, and it is 
clearly the only viable approach when dealing with large 
datasets over wide areas. 

ESIS covers ~ 5 [deg^] in 5 optical bands and is 
based on imaging with the Wide Field Imager (WFI, 
at the focal plane of the 2.2m La Silla ESO-MPI tele- 
scope, iBaade et al.1 Il999l) and the Visible Multi Object 
Svectroqravh (VIMOS . on VLT, iLe Fevre et all 120021: 
iD'Odorico ei~ani2003l) to - 25 - 26 mag in BVRIz. The 
total amount of scheduled observing time is 27 nights with 
the WFI and 8 nights with VIMOS. 

Prime motivations for ESIS are to: 

— obtain optical identification for the roughly 300000 IR 
sources detected by SWIRE in the 5 [deg^] of the ESIS 

^ SWIRE web page: |http : / /swire . ipac. calt ech . edu/ 1 



Berta S., et al.: ESIS 



3 



RA (hours) 

0.8 0.7 0.6 0.5 0.4 



-42 



-43 



-45 



- ' 


1 1 


1 


1 

A 














1 1 


1 


- 






/ 








) ^ 








- 








• 1 


r. 


1^ 




» 1! 


• If 




/ 




- 








m/i 


• 


2 






'*5 






/ 












7' 


• 


I 




"y 


• B 




i 
















• 


) 




i8 


• 7 




















• 


! 




If 


























n 
























WII 











12 10 8 6 

RA (degree) 

Fig. 1. ESIS WFI (BVR) pointings scheduled in the 
ELAIS SI region: a total of 22 fields (red numbered small 
squares) of 30' x 30' are being observed as part of the 
ESIS survey. The region obser ved by ISO llOliver et al.1 
l2000|l . the Spitzer/SWIRE fie ld llLonsdale et alJ2003j) and 
the ATCA radio survey by lGruppio'nret''an lll999') are 
marked. The four circles represent XMM observations. 
Regions 3-6 have also been observed with SOFI/NTT with 
J and Ks filters, during Summer 2002. See Section|H|for a 
review of the available data in ELAIS-Sl. 

area; the present statistics on Spitzer sources indicate 
that ~ 80% of the SWIRE IRAC sources could be 
detected to S = 26 and V = 25.5; 

— provide colors and rough morphologies for source clas- 
sification; 

— build the basis for photometric redshifts for all IR 
sources, and optimize the spectroscopic IR and optical 
follow-up; 

— provide UV-blue restframe luminosities of IR galaxies; 

— produce independent optical samples, selected on the 
basis of their restframe UV-blue properties, for com- 
parison with the IR-selected ones. 

This paper describes WFI observations, data reduction 
and analysis in the central ELAIS-Sl 1.5 [deg^]. VIMOS 
observations will be characterized in a forthcoming paper 
and subsequent releases of WFI data will be presented on 
the ESIS official web page. The data described here are 
being included in the third SWIRE Spitzer Legacy data 
release (Fall 2005). 

Section 2 deals with optical Wide Field Imager ser- 
vice mode observations. Data reduction is described in 
Section 3, where astrometric accuracy, catalog extraction 
and photometric calibration are also discussed. In Section 
4, we present a set of simulations carried out with the 
purpose of estimating the effective depth of the survey 
and flux uncertainties. Section 5 deals with source num- 
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Fig. 2. Transmission curves for ESIS WFI observations. 



Exp. # 




A5 


1 








2 


+2,2 


+24 


3 


-32 


-h48 


4 


-48 


-24 


5 


-h48 


-48 



Table 1. Dithering pattern in [arcsec] of WFI BVR obser- 
vations. Each pointing consists of 30 dithered exposures 
(6 sets of 5). The shifts reported here refer to the first 
exposure of each set. 



ber counts, disentangling point-like and extended objects 
and comparing ESIS to literature data and predictions of 
galaxy evolutionary models. In Section 6, the ESIS cata- 
log is matched to multiwavelength data in the ELAIS-Sl 
area, spanning from the X-rays to the far-infrared, illus- 
trating the potential of panchromatic studies of galaxies. 
Finally, in Section 7, we summarize current findings. 

2. Observations 

The ESIS WFI survey consists of BVR imaging observa- 
tions of ELAIS-Sl, carried out in Service mode between 
October 2001 and January 2005. Further observations are 
scheduled and will be continued until completion of the 
entire 5 [deg^] planned area. 

The ESIS field has been divided into 22 different sub- 
areas, each corresponding to one WFI 30' x 30' pointing. 
The different pointings overlap by about 1.5' on each side 
to obtain a full coverage and to allow trimming individual 
images to avoid edge effects. 

Figure H shows the WFI coverage of the ELAIS-Sl 
ESIS sky area, as well as the SWIRE, ISO, XMM, ATCA 
fields. The SWIRE and ESIS fields are shifted relative to 
one another because of resizing of the original Spitzer area. 
Section El will describe the available multi- wavelength ob- 
servations in ELAIS-Sl. 

The WFI array is made of 8 CCDs, separated by small 
gaps. In order to fill these gaps and correct bad pixels and 
cosmic ray events, each 30' x 30' field is covered by 6 sets 
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99 
1.00 
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9s /in /n9 
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27/10/02 
30/11/03 
28/11/03 
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1.120 
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1.418 
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0.70 


21/10/03 


1.141 


0.78 



Table 2. Summary of WFI, BVR observations. Each set of exposures consists of five 300 [s] frames. Observing dates 
and average airmasses are reported for each set. Seeing values have been measured directly on science frames. 



(Observing Blocks, OBs) of 5 dithered exposures. At least 
one of the six OBs per pointing per band was observed 
during a photometric night, as required by the schedul- 
ing. The standard WFI dithering pattern (see Tab.QJ was 
used. The scale of the instrument is 0.238 [arcsec/pixel]. 

Each single exposure consists of 300 [s], for a total 
exposure time of 2.5 hours per pointing per band (B, V, 
R). 

The filters used are WFI B/99, V/89 and Rc/162. 
Observing Blocks defined after September 2002 (i.e. start- 
ing from pointing no. 5) have B/123 instead of B/99, 
because the standard WFI filters and calibration plan 
changed; in this paper we refer to B,V,R, regardless of 
which are involved. There are no significant differences 
in the final catalog, which is calibrated to the Johnson- 
Cousins system (see Sect. 13. 4f) . Figure 12 shows the trans- 
mission functions for the 4 filters. 

Table 121 summarizes the observing logs for the central 
6 fields (i.e. 1.5 [deg^]). Field 6 observations have not been 
completed yet, but this field is included in this paper be- 
cause it covers part of the area observed in the X-rays and 
near-IR. The seeing and airmass values are the averages on 
the 5 images constituting each OB. The majority of nights 
satisfied the requested observing constraints: seeing<1.2", 
moon distance >120°, moon phase < 0.4. Aborted, or out- 
of-requirements OBs are not listed in Table |21 apart from 
those for pointing no. 6. 



3. Data Reduction 

The reduction of ESIS data was performed within the 
IRAF^ environment, using the package WFPRED devel- 
oped by two of us (LR, EVH) at the Padova Astronomical 
Observatory. This package consists of a set of pro cedures 
that access and upgrade the standard MSCRED I Valde j 
Il998f l tasks. 

Correction of bias and flat-fielding followed the stan- 
dard procedure. Sky flat-field frames were acquired during 
each night and applied to images obtained on the same 
date. 



3.1. Super-Sky-Flat Field 

Even if sky flat-field frames were used, significant illumi- 
nation gradients exist on the pre-reduced frames. Hence 
a further correction was needed, both to obtain a uni- 
form background and to build frames with the same pho- 
tometric zeropoint across the whole field of view. The main 
causes for these gradients are: 

1. sky flat fields point at different positions on the sky 
than the science target, and therefore have different 
relative positions to the moon and very bright stars; 



The package TRAP is distributed by the National Optical 
Astronomy Observatory which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Fig. 3. Example of R-band super-sky-flat frame (negative 
rendered). 




Fig. 4. Weight mask for R-band ESIS WFI pointing no. 
2, used during the detection of sources and the estimation 
of photometric uncertainties. 



2. bright off-axis sources produce spurious multiple re- 
flections in the focal reducer of WFI; 

3. the moon, when present, produces smooth gradients in 
the background of science frames. 

To address these issues, Super-sky-flat field frames (here- 
after ss-flat) were produced, by combining all science 
frames obtained during a given night. The adopted proce- 
dure is sketched below. 

— In order to create the ss-flat, astronomical objects 
must be removed from the science images. Detection 
of all sources and build up of object masks were per- 
formed with the obj masks task in the MSCRED IRAF 
package. 

— Very bright sources (e.g. luminous stars and big galax- 
ies) produce very extended halos on the images. All 
science frames were visually inspected and all halos 
were manually masked. 

— All science frames were then combined together, us- 
ing object-masks to cancel-out the light coming from 
astronomical sources. 

— The result was fitted with a 2-dimcnsional 4th order 
Legendrc function, obtaining a smoothed, but accu- 
rate representation of the ss-flat. Alternatively, a sim- 
ple boxcar smoothing was also tested, but it was found 
that Legendre fitting led to better results, because it 
avoids "holes" produced by the large masks on ex- 
tended halos. 

An example of R band ss-flat is shown in Figure |21 

Since the targeted fields lie all roughly at the same 
equatorial coordinates, the same ss-flat was applied to all 
the images taken during a given night, unless significant 
changes in sky conditions (e.g. presence of atmospheric 
cirrus) occurred. 



3.2. Astrometric calibration and co-addition 

Reduction of wide field imaging data requires a detailed 
astrometric calibration, in order to take into account the 
intrinsic distortions of the instrumentation, and minimize 
the effect of projecting a large sky area (intrinsically spher- 
ical) onto a planar CCD. 

In the case of WFI ESIS data, the 30' x 30' field of view 
requires particular care when mapping pixels into celestial 
coordinates. In view of multi-wavelength cross-correlation 
and future spectroscopic follow up in the ELAIS-Sl area, 
as good an astrometric calibration as possible is necessary. 

An accurate solution is provided by the TNX transfor- 
mation, that combines a linear projection of the sky sphere 
onto a tangent plane (the standa rd gnomonic algorithm, 
TAN, ICalabretta fc Greisenll2002^ and a polynomial func- 
tion for distortions. A simpZi/ied description of the adopted 
mapping is: 



f lngcor[£„T]] 
\ latcor[£_, rj] 



f ail 


ai2 \ 


r 




V 0'2l 


a22 J 




-Vc ) 



(1) 

(2) 



where (^', rj') are the celestial coordinates on the sky, 
mapped by the TNX algorithm, (^, 77) represent the first- 
order solution of the TAN transformation (given by the a 
matrix), {x,y) and {xcUc) are the cartesian coordinates 
of the generic pixel and of the image center, Ingcor and 
latcor are the polynomial equations mapping distortions 
along longitude and latitude. 

The solution of the a bove equatio n array is built on 
dedicated observations of IStone et alJ (|^99) astrometric 
fields D and E, carried out as part of the ESIS program. 
Since observations are spread along a wide period of time, 
several different astrometric solutions have been necessar- 
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Fig. 5. Final, co-added, R-band frame of ESIS WFI Field number 3. 



ily built, because extra-ordinary maintenance of the WFI 
camera affected the pixels— >■ celestial coordinates map. 

The best results were obtained using 4th order distor- 
tion polynomials. Every individual image was astrometri- 
cally calibrated before co-addition. After the solution was 
applied, frames were finely re-centered by ni eans of rigid 
shifts to the GSC 2.2 (STScI & Oa,TO. l2nmh catalog. 

During co-addition, all frames were re-projected, 
transforming their astrometric map to a common one 
(that for the B band), regardless of the photometric band. 
This choice allows the B, V, and R images to be eas- 
ily compared and 3-band catalogs to be straightforwardly 
produced. 

Each one of the 22 ESIS WFI pointings consists 
of at least ~30 science frames per band (B,V,R, see 
Table OJ. While co-adding all these frames together, we 
have excluded all those images with seeing larger than 
~ 1.2 [arcsec], with the exception of pointing no. 6, for 
which few images are available. The FWHM of point-like 



sources (see next Sections), as measured by SExtractor 
(Bertin & Arnouts 199^ on the final images, is reported 
in the fourth column of Tab. 13 

Figure 01 shows the weight map obtained during co- 
addition of Field no. 2, R band. The darkest regions have 
about half the exposure coverage of the lightest ones, due 
to gaps between the CCDs and dithering. The weight value 
defined as the ratio Nimgs/Ntot between the effective num- 
ber of images covering the given pixel and the total num- 
ber of frames belonging to the given ESIS pointing. The 
weight maps were used by SExtractor in object detection 
and noise calculation. Figure [S] shows the R-band image 
of WFI Field number 3; Figure El includes some zooms on 
colorful sources. 

The coordinates of the objects detected on the final 
co-added R-band image of one WFI field are compared 
to GSC 2.2 sources in Fig. (left panels). The r.m.s of 
the distribution of coordinate differences turns out to be 
~ 0.12 [arcsec] in a ■ cos (5 and 0.10 [arcsec] in 6, similar 
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to wha t found bv lArnouts et alJ l)200l[) and lMomanv et all 
l)200l[) . The top left panel of Fig. [7| includes also com- 
parison of data to Gaussian distributions with the same 
variance a. The central plot shows that 70, 80 and 90% 
of the sources are included within A of ^0.13, 0.16 and 
0.20 [arcsec] respectively^. Similar results are obtained for 
the B and V bands and the other pointings. In the bottom 
graph, coordinate differences are plotted against (a, 5): no 
systematic trends are detected. 

Relative astrometry (i.e. the difference between coor- 
dinates in the 3 bands) is very accurate, because the three 
monochromatic frames were registered to the same astro- 
metric map, as mentioned above. The right panels of Fig. 
[3 show the behavior of B vs. R coordinates for one WFI 
pointing; r.m.s are below 0.05 [arcsec]. Similar results are 
obtained also for V band and the other fields. 

The astrometric accuracy has been checked also in the 
regions of overlap between ESIS fields. The bottom right 
panel of Fig. [3 reports the result for GSC 2.2 sources: 
coordinates measured on each pointing are compared to 
those estimated on contiguous fields. The plot includes all 
the central six ESIS fields. 

3.3. Catalog extraction 

Source extraction and magnitude meas urement were per- 
forme d using the SExtractor software ijBertin fc Arnoutl 
Il996l) . 

Thanks to the astrometric transformations described 
in Section the B, V and R images of each WFI point- 
ing are perfectly aligned (see Fig. EJ) with each other. 
Therefore it is very easy to cross-correlate single band 
source lists and generate a 3-band catalog. 

We ran SExtractor on a B-fV-|-R image obtained by 
simply summing the three monochromatic frames belong- 
ing to each WFI field. A Gaussian-filtered image was used 
for detection, with a kernel matching the seeing. We use 
this first 3(7 detection on the BVR image to build a list of 
objects (ASSOC_LIST) to be extracted and measured on 
the individual monochromatic frames. In this way, spuri- 
ous detections on individual images are also minimized. 

Photometry was then performed at 3f7 above back- 
ground r.m.s. with SExtractor also. As a consequence 
of the dithering observing technique, regions along image 
borders and CCD gaps are covered by a small number of 
frames; weight maps (see Sect. l3.2|l are used by SExtractor 
while detecting sources and computing magnitude uncer- 
tainties. Nonetheless, it is worth noting that magnitude 
uncertainties (<0.1 mag down to magnitude ~25, for all 
the three bands) computed by SExtractor likely underes- 
timate the real errors in photometric measurement, be- 
cause they include only photon noise statistics (see tests 
in Section 

Table |31 reports the number of sources detected and 
measured by SExtractor on ESIS individual WFI frames. 
The total number of sources in pointings 1 to 6, is 132712, 

3 defined as A = cos(5)]2 + [A(5)]2 
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Table 3. Summary of source extraction. For each of the 
6 ESIS WFI pointings analysed here, the total number 
of co-added frames, number of 3tT extracted sources and 
FWHM of point-like objects are reported. 



smaller than the sum on individual images because of over- 
lap; the total number of objects on Fields 1-5 — full depth 
— is 118197. 

3.4. Photometric calibration 

The fundamental requirement to be fulfilled by the ESIS 
service mode WFI program consisted of observing at least 
one OB (i.e. one set of 5 exposures, hereafter ^^refer- 
ence OB") per pointing, per band, during one photomet- 
ric night. During the same photometric night, imaging of 
Landolt (1992) Rul49 or TPHE spectrophotometric stan- 
dard fields — at similar airmasses — was performed. Eight 
different exposures were taken each time, in order to in- 
clude the main standard stars on each individual CCD of 
the WFI array. 

When combining ~30 different images, belonging to 
various nights (see Tables |21 and |3J) , frames obtained with 
different sky conditions are mixed, hence the photomet- 
ric information gets lost. Having a photometric reference 
night is necessary, to recover a correct calibration of mag- 
nitudes. The five images belonging to the reference OB 
were combined together and were used to determine the 
zeropoint shift Amag=magcoadd— magre/, caused by the 
co-addition of ^30 frames into the final frame. 

The preliminary catalog fSect. lB^ was matched to the 
photometric conditions and normalized to unit airmass 
and exposure time, using the equation: 

magp/jot = mag 

coadd + 2.51ogio(tea:p) (3) 

-Kx ■ A.M. - Amag 
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Fig. 7. Astrometric accuracy of ESIS WFI data. Left panels: absolute R-band astrometric accuracy. The equatorial 
coordinates of ESIS sources are compared to those in the GSC 2.2 catalog. The central plot shows the distribution 
of a/[A(q; • cos(5)]2 + [A((5)]2; the vertical hues include 70, 80 and 90% of the sources respectively. In the bottom left 
panel, the values of A(a • cos 5) and A{S) are plotted against {a, 5). Right panels: relative astrometric accuracy. The 
coordinates measured on the B band image are compared to those in the R band. In the bottom right plot, the difference 
in coordinates among contiguous fields is shown, for common GSC 2.2 sources. 
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Fig. 8. Self-check of photometric caUbration. The color 
equations (Eqs.0J built on three photometric bands, lead 
to two different solutions, based on the {B—V) and {V—R) 
colors. The magnitudes based on the two calibrations turn 
out to be perfectly consistent to each other. 
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where A.M. is the average airmass of the reference OB, 
and K\ are the atmospheric extinction coefficients pro- 
vided by ESO (ifbgg = 0.23, Kbi23 = 0.22, i^„89 = 0.11, 

Krl62 = 0.07). 

After pre-reduction and ss-flat correction of standard 
fields, color equations to transform instrumental magni- 
tudes into the Johnson-Cousins system were computed for 
each CCD. The results are consistent with ESO^ official 
equations: 

Bj = 699 + 24.65 -I- 0.24 x {B., - Vj) 
Bj = 6123 + 24.71 + 0.19 X {Bj - Vj) 

Vj = VS9 + 24.15 - 0.13 X {Vj - Rc) (4) 

Vj = VS9 + 24.15 - 0.09 x {Bj - Vj) 
Rc = rie2 + 24.47. 

The 699 or 6123 equation is adopted depending on which 
filter was used (see Sect.|21l. The catalog was calibrated to 
the standard Johnson-Cousins photometric system {Bj, 
Vj, Rc), by solving the above array iteratively. 

Since 3 bands and 2 colors are available, two different 
calibrations are possible, based on [B — V) or (V — R). 
Figure |H1 shows the difference in magnitude between the 
two different calibrations. When one color was not avail- 
able, a spiral-like average of 1.0 was adopted. Only B and 
V are shown in Fig. |S1 since no color term affects the R 
band filter. The two estimates are self-consistent; the as- 
sociated magnitude uncertainty is similar or even smaller 
than that estimated by SExtractor on background noise. 

see ESO-WFI homepage 



Fig. 9. Detection efficiency of synthetic sources in the V- 
band. Three cases are studied: a real science frame already 
containing real objects, a ss-flat image and a totally syn- 
thetic sky frame (see text for further details). Histograms 
represent the real data, the continuous red line is a spline 
interpolation. Vertical dotted lines set the 97% and 90% 
limits. The bottom panel includes the difference in detec- 
tion rate between tlie three cases. 

4. Quality tests 

In order to test data quality and survey performance, a 
series of tests have been performed, based on simulations, 
aimed at determining: 

— a semiempirical estimate of the effective depth reached 
and of completeness; 

— a reliable estimate of magnitude uncertainty, overrid- 
ing SExtractor's; 

— the reliability of SExtractor's stellarity flag. 

To this end, we have added synthetic sources to our im- 
ages, by using IRAF tasks in the package artdata. 

The analysis comprises two steps. First we have tested 
three different types of images in order to understand the 
best way to reach our goal: 

i) the real image, containing all real sources; we chose 
WFI pointing number 2, because it represents the 
worst case, containing a very bright (V ~ 7.8) star; 

ii) an empty image obtained from the non-smoothed ss- 
flat frame, built from 30 individual images belonging 
to the 6 ESIS central pointings; 
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iii) a totally synthetic background image generated by 
IRAF, matching the r.m.s. properties of the real im- 
age. 

This test was run on the V band and using point-like ob- 
jects only; results are reported in Section 14.11 The best 
choice turned out to be the ss-flat frame, because it re- 
produces all the defects of real science images, but is not 
affected by confusion problems. 

A number of different simulated images per band were 
then produced, by adding point-like, De Vaucouleurs or 
exponential disk objects onto the ss-flat. In each case, a 
different image was produced per 0.25 magnitude bin, in 
the range 20-27 mag. A population of 2000 sources was 
added in each case, for a total of '--^50000 objects on 14 
different images per band. 

Concerning De Vaucouleurs and disk profiles, intrin- 
sic half-light radii between 2 and 10 [kpc] and a Euclidean 
power law luminosity function were adopted. The artdata 
package assumes a redshift distribution that corresponds 
to the apparent magnitude distribution. The standard cos- 
mological dimming of flux and angular size are applied 
to each artificial galaxy. Finally the synthetic objects are 
convolved with a gaussian kernel consistent with the real 
seeing of observations and are added to the ss-flat frames. 

The results described below are based on SExtractor's 
performance on this simulated images. 

4.1. Detection rate 

Figure |51 summarizes the results of the preliminary test 
performed on the V band images. The detection efficiency 
of SExtractor was defined as e = NsEx/Nin, i.e. the frac- 
tion of input synthetic sources recovered by SExtractor. 
The three top panels in Fig. IHlshow the dependence of e 
on magnitude, measured on the real science final image, 
on the ss-flat frame, and on a synthetic sky image. The 
two bottom panels report the difference in e, between the 
various cases. 

The effects of confusion influence the detection of syn- 
thetic sources on the real science frame, with respect to 
other cases, in two ways: e decreases from unity at brighter 
magnitudes, because sources get lost in the halos of bright 
stars and very extended galaxies; e decreases more slightly 
than in the other cases, because of spurious detection due 
to superimposition onto real sources. 

Concerning the ss-flat image, which was intentionally 
neither smoothed, nor cleaned of bright halos, e decreases 
less suddenly than on the synthetic sky, because some spu- 
rious sources are detected in the halos. 

We estimate the amount of spurious detections by run- 
ning SExtractor on the empty ss-flat (i.e. without any 
simulated object added), searching for sources at the po- 
sitions of those previously detected. No spurious sources 
were detected down to V~25; at fainter fluxes about 1-2% 
of sources are not real, in the range 25-26 mag. Similar 
results were obtained in B and R. 
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Fig. 10. Top panel: V-band ss-flat simulation. Difference 
between the input magnitude of point-like synthetic 
sources and the value measured by SExtractor, as a 
function of input mag. The solid line traces the me- 
dian, the short-dashed standard deviation, and the long- 
dashed semi-inter-quartile (s.i.q.). Bottom panel: compar- 
ison between the r.m.s. (or s.i.q.) computed on simulated 
data and the median magnitude uncertainty measured 
by SExtractor on real data, in a given magnitude bin. 
Solid, dotted and dashed lines represent point-like, disk 
and spheroidal sources, respectively. 



Table 01 summarizes the resulting completeness values 
derived for point-like, De Vaucouleurs and exponential- 
disk simulations on the ss-flat images, in the three 
bands. Overall, a 90% detection efficiency is reached at 
mag'^25.5, with some scatter, depending on the case. For 
pointing number 6, consisting of 10 exposures only, this 
limit is ~ 0.5 mag brighter. 

4.2. Magnitude uncertainty 

An estimate of magnitude uncertainties is obtained by 
comparing the input simulated fiuxes and the measured 
values. The top panel in Fig. 1101 shows the case of point- 
like sources added to the V band ss-flat frame: the differ- 
ence between input and extracted magnitudes is plotted 



12 



Berta S., et al.: ESIS 



Point-like sources 





R 


V 


B 


mag 




s.l.q. 


rate 


(7 


s.l.q. 


rate 


a 


s.l.q. 


rate 


21 


0.024 


0.005 


0.999 


0.024 


0.003 


0.999 


0.022 


0.002 


0.999 


23 


0.063 


0.034 


0.999 


0.042 


0.020 


0.999 


0.033 


0.015 


0.999 


25 


0.283 


0.175 


0.999 


0.199 


0.115 


0.999 


0.158 


0.095 


0.999 


26 


0.478 


0.266 


0.064 


0.291 


0.178 


0.953 


0.282 


0.172 


0.980 



Exponential-disk sources 
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De Vaucouleurs sources 
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Table 4. Summary of quality simulations on the ss-flat image. Three cases are reported, consisting of point-like, De 
Vaucouleurs and exponential-disk synthetic objects. For each band, magnitude standard deviation and semi-inter- 
quartile (s.i.q.) are listed, as well as extraction efficiency. 
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Fig. 11. Left, dependence of the stellarity index on Vj magnitude, for simulated point-like sources, simulated expo- 
nential disks and real data. The horizontal dotted line sets the CLASS_STAR=0.95 limit. Right: half-flux radius as a 
function of Vj magnitude. The dark dots define the locus of point-like sources. 



as a function of input magnitude. The solid line represents 
the median difference, the short-dashed line is the stan- 
dard deviation from the median value and the long-dashed 
line traces the semi-inter-quartile^. 

The bottom panel in Fig. EH compares the standard 
deviation and semi-inter-quartile of simulations (see top 
panel) to the median magnitude uncertainties measured 
by SExtractor on real data. The latter is computed as the 
median of uncertainties for all real sources in a given mag- 
nitude bin. In the case of point-like sources, SExtractor's 
errors are compared both to the standard deviation from 
the median (filled circled) and to the semi-inter-quartile 

^ defined as (53 — qi)/2, where qi and qs are the first and 
third quartiles. The first quartile is the number below which 
25% of the data are found and the third quartile is the value 
above which 25% of the data are found. 



(open circles). The synthetic a is typically driven by out- 
liers (such as unresolved double stars), while the semi- 
inter-quartile width of the distribution provides a fairer 
comparison to SExtractor's uncertainties. 

The estimate of magnitude errors for non-pointlike 
objects is compared to the simulated ones for disk and 
spheroidal galaxies (open squares and triangles respec- 
tively). The discrepancy between the two estimates is 
larger in this case, because of non-shot-noise effects like 
failure in detecting the faint wings of galaxy profiles or 
size dimming for the faintest sources. 

Table^lsummarizes the results for all the cases studied. 
Typical uncertainties at mag= 21, 25 are ^ 0.03, 0.3 in 
the B band, ~ 0.05, 0.5 in V and R. An uncertainty of 
0.05 mag at 21 corresponds to one of 5% in flux, while 0.5 
at 25 corresponds to 60% in flux. 
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4.3. Stellarity 

Finally, the simulations described above are also useful to 
test the reliability of SExtractor's stellarity flag. In com- 
bination with real data, it is possible to evaluate down 
to which magnitude the CLASS_STAR flag represents a 
realistic indication of the source profile. 

The left plots in Fig. ^2 show the distribution of 
CLASS_STAR as function of magnitude, for the simulated 
point-like sources (left panel), simulated exponential disks 
(central panel), and for the real data (right panel). 

The plots suggest that a CLASS_STAR>0.95 is a good 
representation of point- like sources up to 23. Note 
that effectively no point-like sources are detected above 
this threshold, at fainter fluxes. Nevertheless, a significant 
contamination by galaxies is expected at fluxes fainter 
that V~21, as confirmed by Fig. El on real data (see also 
Section l5. 1(1 . Similar results are obtained in the other two 
bands, with slightly fainter magnitude limits. 

An intrinsically wide scatter in CLASS_STAR is typi- 
cal of galaxies. Moreover these tend to behave like point- 
like sources at faint fluxes, because the profiles of small 
objects are dominated by seeing. 

A safer way of extracting point-like sources relies on 
the half-light radius, ri/2, enclosing 50% of the object's 
total flux. As for point-like sources ri/2 depends only on 
the image seeing, when plotted against magnitude it de- 
fines a tight "stellar" locus (see the diagram on the right 
side of Fig. Ill() . Although still contaminated by QSOs, 
this method provides a good representation of point-like 
objects down to 22. 



ri/2 provides a good representation of point- like sources, 
but this method is limited to magnitudes brighter than 
~ 22. Moreover, it is worth to note that quasars are in- 
trinsically point-like sources and cannot be a priori distin- 
guished from stars, on the basis of BVR colors only. 

The long-dashed line in Figure El represents the ex- 
pected magnitude distribution for s tars, based on the 
Milky Way model by IJarrettI l)l992|) and Ijarrett et all 
|199J). These authors adopt a three component model, 
including a halo, a thin disk and a thick disk stellar pop- 
ulation. This model reproduces the observed star counts 
in different galactic lines of sight, as well as SDSS data. 

The Milky Way model is consistent with ESIS data 
at bright fluxes, in the B and V bands, while it seems 
to underpredict (by about 20-25%) star counts in the R 
band at magnitudes < 19. Between magnitudes 21 and 
22, observed pointlike counts still show some excess with 
respect to Jarrctt model. This may be caused by resid- 
ual contamination from faint unresolved galaxies or by a 
quasar component. 

In what follows (e.g. Figure [TS|l . two corrections will 
be applied to numb er counts for stars: one based on the 
I.Tarrett et al.l l(l994l) model, the other using SExtractor's 
identification. Overall the latter correction is more solid 
than the former at bright fluxes, while the former is most 
suitable for faint magnitudes. In any case it should be 
noted that the stellar contribution to the counts at the 
faint end (less than 3%) is almost negligible, when com- 
paring data and galaxy evolutionary models, while it dom- 
inates the bright counts. 



5. Source counts 

The magnitude distribution of galaxies, per unit sky area, 
provides strong constraints on th eir evolutionary history 
l Ellijll997t iMetcalfe et al lll995() . In monolithic-collapse 
models ((Eggen et all 11961 , the bulk of stars formed in 
early- type galaxies at very high redshift, and galaxies 
evolve passively to present day; while subsequent merging 
and star formation are limited. In a h ierarchical merger 
scenario l|Kauffmann fc Chariot! 1 199 j^ . the bulk of the 
stars form in disk-like galaxies that later merge to be- 
come early- type galaxies. The luminosity and density of 
galaxies strongly evolve with cosmic time due to mergers 
and related starbursting events. 

We build ESIS counts excluding WFI pointing number 

6, hence including sources within a 1.22 [deg^] area. Figure 
ll2lshows the data: the solid thick histogram represents the 
total observed number counts, the dotted histogram the 
counts for point-like sources (see Sect. 14.31 . and the thin 
solid histogram the counts for extended objects. 



5.1. Star-galaxy separation 

In order to compare to galaxy evolutionary models, stars 
must be removed from number counts. On the basis of 
the analysis presented in Sect. 14.31 the half-light radius 



5.2. Comparison to literature data 

Figure El also compares ESIS galaxy counts to a 
collection of l iterature data. In particula r, HDFN 
llWilliams et a,lJ (l99(i) . HDFS ( Metcalfe et a,].l 12001 h . 
SPSS (lYasu da et a,lJ EoOli) . VIRMO S deep survey 
((McCra cken et al.l l2003l)^^aCD F (Alc al^_^et all l2004l) . 
NTT SUSI D eep Field lArnouts et al. 1999). Additional 



data are from IMetcalfe et all lll99lll200L) . .Picardl l)l99l|) . 
Wilson (2003) and Gardner et alJ l|l9!m All magnitudes 
have been transformed to the Johnson- Cousins photomet- 
ric system, adopting the appropriate relations (usually re- 
ported by authors). 

ESIS number counts are generally consistent with the 
results from other surveys, in the common flux range. 
At bright magnitudes, some differences arise between our 
counts for extended sources (which can be considered 
galaxies for mag< 18). Th is is particularly true in the V 
band, the Gardne r et alJ (^996) data being significantly 
higher than ESIS', possibly due to a different choice in 
accounting for stars, or to cosmic variance. 

5.3. Galaxy evolutionary models 

Differential galaxy number counts can be described as the 
integral over cosmic time of the galaxy luminosity func- 
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Fig. 12. Observed ESIS number counts (histograms), compared to literature data (dots) and Ijarrett et all l|l994) 
Milky Way model (long-dashed line). The ESIS data are split into extended (this solid histograms) and point-like 
(dotted) components. The bold solid histograms represents the total ESIS counts. 



tion (LF), in a specific flux bin (e.g. ^anceschini et alJ 
l200lHGuiderdoni fc Rocca-Volmerangelll99l[l . The shape 
and normalization of the LF is assumed constant in non- 
evolutionary models; galaxy luminosity evolves by means 
of galaxy ages, in pure luminosity evolution models (PLE); 
finally in complete evolutionary models, the LF changes as 
a function of redshift, for example due to enhanced star 
formation activity, triggered by encounters and mergers 
between galaxies. 

In the recent years hierarchical clustering models have 
gained increasing popularity in the description of galaxy 
formation and evolution. Primordial density perturba- 
tions are amplified by gravitation and collapse to form al- 
most virialized structures called dark matter (DM) haloes. 
In the cold dark matte r (CDM) scenario (Peebles .198 j; 
iBhimenthal et alJll98^ . smaller haloes form first, subse- 
quently merging into bigger haloes. Gas cools down in 
the potential wells of the haloes, leading to the forma- 
tion of stars and galaxies. The spectrophotometric evolu- 
tion of stellar populations, combined with the history of 
stellar formation of galaxies, finally produce their observ- 
able properties: spectra, colors, etc. Semi-analytic modern 
models combine these and others ingredients and follow 
the various processes transforming the primordial power 
spectrum of density fluctuations into the spectral en- 
ergy distributions of stellar populations and galaxies, (e.g. 
[Hatton et al. 2003; Guidcrdoni et al. 1998). 

Figure El compares ESIS galaxy counts to the four 
different kinds of models. Galaxy counts are obtained 
from total number counts in two ways: by subtract- 
ing point-like identifica, tions (open squares) and by using 
the lJarrett et alJ l)l994|) model (stars). Poisson errors are 
smaller than the symbols, for mag>18. 

Completeness correction is attempted at the faintest 
fluxes, on the basis of the detection rate analysis carried 
out in Section 01 Only one magnitude bin can be fully 
recovered, reaching B,V~24.5 and 24. Table |S1 lists 
the ESIS counts. 
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Table 5. ESIS BVR number counts [deg"^ (0.5 mag)"!]. 
Completeness corrections has been applied, based on sim- 
ulations of galaxy-like sources (see Sect.0J. Total and ex- 
tended source counts are provided. 



Describing the data as a power law 

log N (X -f X mag, (5) 

we find = 0. 48 in the range B = 1 9 — 24, similarly to 
lEIlij (|1997^ a,nd iMetcalfe et a1.l ll2nnih . The V-band data 
can be fitted with a double power law, having an elbow at 
V - 20: 7v = 0.56 between 16 and 20, and 7y = 0.37 for 
= 20 — 24. Finally, concerning R band data, the galaxy 
counts are well reproduced by 77? = 0.58 in the range 
16 - 20 mag and 7r = 32 for i? = 20 - 24, consistent 
with lMetcalfe et alJ l|200l|) . Models are compared to data 
only for B and R band, because — in the literature — the 
R band is usually preferred to the V band. 
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Fig. 13. Comparison between ESIS galaxy counts and prediction of evolutionary scenarios. Open squares are number 
counts for extended sources (based on the half- flux radius analysis). Asterisks represent total galaxy counts corrected 
for stars contribution with the J arret t et al. (199^) Milky Way r nodel. Filled squares are d ata corrected for completeness 
(at t he faintest fluxes). M odels belong to lMetcalfe et al.l l|200lj) . lNaEfashima et al.l l 2002l) and the GallCS team (shaded 
area, lHa,tton et ablEnni . 



In Fig. ^1 the long-dashed lin e repr esents the no- 
evolution model, bv iMetcalfe et ahl l)200lh . obtained for 
a go = 0.5 cosmology. These authors model the galaxy 
luminosity function (LF) as a Schechter function, with 
a = —0.7 for E/SO/Sab galaxies, a = —1.1 for Sbc spi- 
rals and a — —1.5 for bluer Scd/Sdm ones. In the B-band 
the model was norrnalized to B ^ 18, instead of 15 (see 
IMetcalfe etallllQQll llQQSl l200li for a discussion in sup- 
port of this choice) and reproduces reasonably well the 
observed counts down to _B 21. In the R band the con- 
sistency with data holds only to R'^ 19.5. At fainter fluxes 
the model quickly diverges from the ESIS counts. 

The shor t-dashed line in Fig. ^| is a PLE model 
ijMetcalfe et al..2 001'). built assuming that the evolution of 
the different morphological types is governed by their star 
formation history. Exponentia lly decaying star formation 
rates are adopted jBruzual fc Charlot.1993) , using t = 2.5 
Gyr for E/SO and Sab galaxies and r = 9 Gyr for other 
spirals types. Present day galaxy ages are ~ 12.7 Gyr, im- 
plying a formation redshift of = 9.9. This model is in 
good agreement with ESIS data down to B,R^23, after 
which the accordance fails. 

IMetcalfe et a,1.ll|2nnih introduced a population of dwarf 
ellipticals to increase the model-predicte d counts a t the 
faintest fluxes (down to B,R~ 28, see also lE11iijll997j) and 
match observations. Alternative solutions lead to equally 
good fits, for example merging evolutionary models. 

The solid line in Fig. 1131 represents an evolutionary 
model, with a LF slope of a = —1.75 for blue Scd/Sdm 

J :alaxies at redshift z > 1, steeper than the local a = —1.5 
Metcalfe et"alll200l|) . This model provides a good fit to 
ESIS data over the entire magnitude range considered. 



The GallCS (G alaxies In Cosmological Simulations, 
iHatton et al1l2003|) project provides a hybrid model for 
hierarchical galaxy formation studies, combining large cos- 
mological A^-body simulations with semi-analytic recipes 
to describe the properties of galaxies within dark matter 
haloes. We have retrieved mock catalogs produced with 
the GallCS model, for a total of 10 different fields of 1 
[deg^], in order to take into account the cosmic variance 
of the mock catalogs. 

GallCS simulations assume a ACDM flat Universe 
(a™ = 0.333, r^A = 0.667, h = 0.667, o-g = 0.88). Each 
realization consists of a 100 [Mpc] cube, containing 
256^^ particles of mass 8.272 • 10^ Mq. The spatial resolu- 
tion of the simulation is 29.29 [kpc]; the evolution of the 
Dark Matter density fleld is followed from z = 35.59 to 
z = 0, through 70 different snapshots. 

Dark matter halos containing more than 20 particles 
(i.e. with M > 1.66 • 10^^ Mq) are identified and their 
merging history trees are then computed. Baryons are 
evolved within these halo merging history trees, follow- 
ing a set of semi- analytic prescription s (the Mock Map 
Facility, MoMaF, iBlaizot et all |200 5>) that account for 
(among other effects) heating and cooling of the gas within 
halos, star formation and supernovae feedback on the en- 
vironment, evolution of stellar populations, metal enrich- 
ment, disk instabilities, tidal stripping, and formation of 
spheroids through galaxy mergers. 

The shaded regions in Fig. represent the locus of 
predictions by the GallCS model, for the 10 different sim- 
ulations considered. In the R band, consistency between 
model and observations holds at the brighter fiuxes, but 
the model seems to overpredict galaxy number counts by 
at least ~30% at magnitudes fainter than ~20. In the B 
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(Bj-Rc)vesa (Bj-Rc)' 



Fig. 14. Observed B,j — Rc color distribution. Left: extended sou rces. Overplotted are: a no-evolution (long-dashed 
line. iRoche et al.lll996(l . a PLE (short-dashed. lRoche et^ll99dl and Ga llCS A-CDM (shaded area. iHatton et alJ 
|2003() models. Right: point-like sources, as compared to lJarrettTtal] l)l994j) predictions for spheroid (dashed line) and 
disk (dotted line) stars. 



band, the discrepancy is more significant, model predic- 
tions overpredicting number counts by a factor of ~ 1.5 
at Bj > 23. 

iNaeashima et al.l (|2002|) presented a ACDM scmi- 
analytic model, which produces merger trees of Press- 
Schechter dark halos, in an ri„j = 0.3, J^a = 0.7, h = 0.7 
Universe. This model includes the effects of dynamical re- 
sponse to supernova- and starburst-induced gas removal 
on size and velocity dispersions, which seems to play an 
important role in dwarf g alaxy formation. The do tted line 
in Fig. represents the *Nagashi ma et al ] ("2002?) model, 
which turns out to be fairly consistent with ESIS data in 
all bands, excepted for some excess in the number counts 
at faint fluxes. 

Overall, the comparison of ESIS optical number counts 
to theoretical predictions favors a scenario in which the 
Universe has experienced an epoch with many encounters 
and interactions between galaxies, possibly at z > 1, when 
the mean restframe UV-B luminosity of galaxies was larger 
than locally, because of enhanced star formation activity, 
likely triggered by mergers. 

5.4. Color distributions 

The distribution of observed colors, for both galaxies and 
stars, provides additional constraints to cosmic evolution 
and galactic models. 

The left panel in Fig. El shows the distribution of 
the Bj — Rc color for extended sources, in different Bj 
mag nitude bins. A non -evolutionary model (long-dashed 
line, iRoche et al]ll996|) systematically underpredicts the 



color counts, i n the 22 < B j < 2 4 bins. The PLE model 
(short dashed, iRoche et al.lll996l ) is more consistent with 
the data, but still does not reproduce observations in the 
faintest magnitude bin. Concerning the two brightest mag- 
nitude bins, only few model predictions of the observed- 
frame colors are available in literature. The GallCS A- 
CDM color histogram seems to be too much peaked and 
to predict more sources than observed (see also Fig. ll3|) . in 
particular overestimating the number of faint blue galax- 
ies. 

The right panel in Fig. 1141 compares the observed color 
distribution of point-hke sources (defined as in Sect. I4.3jl 
to the predictions of the IJarrett et"an l|l994|) Milky Way 
model. The dashed line refers to the spheroidal compo- 
nent, the dotted one to the disk population. Despite model 
predictions are fairly consistent to the observed point-like 
number counts (Fig. [T^. the observed color distribution 
shows an excess of red sources at the faint fluxes and might 
caution against a residual contamination from unresolved 
galaxies. 



6. Potential for multi wavelength investigations 

The ELAIS-Sl field benefits from an extensive follow up 
carried out over the whole electromagnetic spectrum from 
X-rays to radio wavelengths (see also Fig. In this 
Section we match ESIS data to catalogs in other spec- 
tral domains, with the purpose of showing some of the 
potential of multi- wavelength analysis of deep surveys. 
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6.1. Observations in ELAIS-Sl and match to ESIS data 

The European Large-Area ISO Survey (ELAIS, 
iRowan-Robinson et al.ll999H01iver et al.l20(¥il) . consisted 
of a deep, wide-angle survey at high galactic latitudes, 
with the Infrared Space Observatory (ISO). The primary 
ELAIS survey was carried ou t over 12 [deg^], at 1 5 
(with the ISOCAM camera, "Cesarsky et al.' 'l996^ and 
90 Atm (with ISOPHOT, Lemke et a l, 1996): additional 
observations were performed in restricted areas at 6.7 
and 175 um, in collaboration with the FIRBACK team 
jPuget et a].llT99 9). The main survey area was divided 
into three fields in the northern hemisphere (Nl, N2, and 
N3) and one field in the southern hemisphere (SI), which 
is the target of the ESIS survey. 

The ELAIS-Sl field, centered at RA = 00'*34™44".4, 
DEC = -43°28'12".0 (J2000.0), covers a sky area of about 
2° X 2° and includes the minimum in Galactic 100 iim 
cirrus emission in the Southern sky ijSchleeel et alJll998() . 

The ELAIS survey in the SI field reached 1.0, 0.7 
and 70 fmJy] depth in the 6.7, 15 and 90 fim bands 
I Rowan-Robinson et all (2004). No 175 /im observation 
was carried out in ELAIS-Sl. The final 15 fim catalog 
contai ns 736 sources down to ~1 mJy, ~20% of which are 
sta rs ifVaccari et al.ll2005l). 

iGruppioni et al.l (^92) performed a radio (1.4 GHz) 
survey with the Australia Telescope Compact Array 
(ATCA), over 4 sq. deg. in the ELAIS-Sl area, produc- 
ing a catalog (already public) of 652 sources, down to a 
minimum flux density of ^0.2 mJy (5a). The current 1.5 
[deg^] ESIS area contains 268 of the ATCA radio sources 
detected in ELAIS-Sl, ~65% of which have an optical 
counterpart. A deeper survey is in progress by B. Boyle 
an d collaborators. 

iLa Franca et al.l ( 2q3) presented spectroscopic and R- 
band data for 406 15 fim sources in the ELAIS-Sl field, 
over the flux density range 0.5 < S'(15/Ltm) < 150 [mJy]. 
The R band data were obtained by an ESO imaging cam- 
paign with the DFOSC instrument mounted on the 1.54 
Danish/ESO telescope at La Silla (Chile), reaching 95% 
completeness level at R^22.5. 

Spectroscopic observations of the optical counterparts 
of the ISOCAM SI sources were carried out at the 
2dF/AAT, ESO Danish 1.5 m, 3.6 m and NTT telescopes, 
adopting instrumental configurations with ~10 A resolu- 
tion and covering the 4000-90 00 A wavelength range on 
average ijLa Franca et alJl2004j) . 

Spitzer/SWIRE observations in ELAIS-Sl cover 
roughly the whole ISO region and a total sky area of 
~7 [deg^], reaching 5a sensitivities of 3.7, 5.3, 48, 37.7 
and 350 u Jy in the 3.6, 4 5, 5.8 , 8.0, 24 /xm channels re- 
spectively l|Lonsdale et al .120041) . IRAC angular resolution 
spans ~2.5 to 3 [arcsec], while for the MIPS 24 /im band 
it is ~6". Data reduction was carried out by the Spitzer 
Science Center and SWIRE team. Absolute photomet- 
ric uncertainty is '^10% both for all bands. SExtractor's 
Kron fluxes are considered for extended sources, while 
aperture photometry is used for point-like objects, cor- 



rected for aperture losses, as recommended by the Spitzer 
Science Center and SWIRE data release documentation 
((Surace et al.l l2004l) . The astrometric accuracy of the 
SWIRE data products is ~0. 5" (after image reconstruc- 
tion and comparison to 2MASS positions). The SWIRE 
band-merged catalog has been cross-correlated with ESIS 
WFI sources, by means of a simple nearest-match algo- 
rithm, adopting a radius of 1"; 71779 SWIRE sources lie 
in the ESIS WFI current area; 53267 matches (i.e. -74%) 
are found using this criterion. 

About 40% of the ELAIS-Sl area was surveyed by 
BeppoSAX, down to a 2-10 keV sensitivity of ~10~^^ [erg 
cm-^s-i] (Alexander et al. 2001). More recent XMM ob- 
servations target the central area of the ISO field (Puccetti 
et al., in preparation). A mosaic of four deep XMM point- 
ings covers ^--^0.6 [deg^] with the EPIC and two MOS cam- 
eras onboard XMM. The average net exposure time is —50 
[ks], in the energy range between 0.5 and 10 [keV]. WFI 
and XMM sources were cross-correlated if lying within a 
radius of 2" (close to the typical la XMM error box); 264 
out of the 479 XMM sources are matched to ESIS objects. 

The XMM ELAIS-Sl area is the target of further 
complementary observations. As part of ESO Large 
Programme 170.A-0143 (P.I. A. Cimatti, Dias et al., 
Buttery et al., in prep.), the ELAIS-Sl central —0.8 [deg^] 
was the target of near-IR J and Ks imaging with SOFI on 
NTT, reaching J<21 and Ks<20 (Vega). 

VIMOS/VLT spectroscopy was carried out in summer 
2004 and 2005, in 2 x 33 hour observing runs (73.A-0446 
and 75.A-0428; PI: F. La Franca). The whole XMM-NIR 
area was the target of low resolution spectroscopy over 
the wavelength range 5000-9500 A. During the first half 
of the survey, 1200 X-ray and Kg < 18.5 galaxies were 
observed; in period 75 the targets consisted of 1000 X- 
ray sources plus 24/im SWIRE galaxies. R-band VLT pre- 
imaging provides additional data down to i? — 25. Further 
spectroscopy of optically bright [R < 21) X-ray and mid- 
IR sources is being performed at the 3.6m/ESO telescope, 
during Fall 2005. These data will be presented in future 
papers by F. La Franca and collaborators. 

The ultraviolet Galaxy Evolution Explorer (GALEX, 
iMartin et alJl2005|) Deep Imaging Survey (DIS) is observ- 
ing —80 fdeg^j in 12 di fferent areas, including ELAIS-Sl 
l)Burgarella et al.ll2005l) . Four ESI GALEX tiles are al- 
ready available to the public as part of data Release^ no.l 
(GRl); two of these (namely tiles 00 and 01) match —0.93 
[deg^] of the ESIS- WFI area described here. The expo- 
sure times for these two tiles are —100 and —50 [ks] re- 
spectively. The photometric catalogs in the far-UV (1344- 
1786 A) and near-UV (1771-2831 A) have been matched to 
ESIS data. In the ESIS-GRl common area, 21025 GALEX 
sources were detected. The in-fiight measured point-source 
FWHM is 4'.'0 and 5'.'6 in the two channels respectively. 
To be conservative, we used a matching radius of 2'.'0 only, 
which is sufficient for our demonstrative purposes of show- 



^ GALEX Release 1 (GRl) DIS data were released on Jan 
4th, 2005 ^ http://www.galex.caltech.edu/ J 
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Fig. 15. Color-color plots for ESIS sources matched to SWIRE/Spitzer objects. Template tracks represent k-corrected 
colors for diff erent classes o f sources, as a function of redshift (z = — 3). The class-V stellar track was obtained by 
extending the lPickle^ ( 199?^ ) library to the infrar ed, assuming a Rayleigh- Jeans law. The reddening arrow refers to the 
standard extinction law by Cardelli et al.l ( Ji 



ing the potential of observations in ELAIS-Sl. In this way 
spurious and multiple matches are avoided. A more accu- 
rate match will be performed in the future, in order to 
provide as complete a multi-A study as possible. In this 
way, 17010 matches were found (~80%). We use the flux 
measurement provided by GRl. 

6.2. Optical-IR colors 

Color-color plots, built on a wide wavelength baseline turn 
out to be very useful to disentangle various classes of 
sources, different spectral domains being sensitive to dif- 
ferent emission components. 

Figure shows three color-color plots including opti- 
cal WFI, and mid-IR IRAC and MIPS photometric data. 



Synthetic tracks are overlaid on the data, as obtained by k- 
correcting template spectral energy distributions (SEDs) 
from 2 = 3 to now. The solid curve belongs to a 5 Gyr 
old elliptical model (truncated at z = 1.2, to be consistent 
with Universe age^), the dotted line is a spiral M51-like 
model, the short-dashed line is a M82 prototype starburst 
templ ate: these th ree models are built with the GRASIL* 
code l|Silva et al ] 11998), and then upgraded by introduc- 
ing observed PAH mid-IR features. The long-dashed line 
belongs to a semi-empirical SED of the blue peculiar spi- 
ral model of NGC 4490 (PoUetta et al., in preparation) 



tH =il3.6 Gyr, for a 0^=0.27, nA =0.73. Ho=71 [km s"^ 

Mpc"^] cosmology jSpergel et al.ll200^^] 

^ GRASIL homepage: http://web.pd. Eistro.it/gr anato/grasil/| 
grasil.html 
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and the dot short-dash model belongs to the type-I AGN 
Mrk 231 (Fritz et al., in prep.). Finally the dot long-dash 
line bel ongs to class V stars and is obtained by extend- 
ing the IPicklesI l)l998l) models to the infrared assuming a 
Rayleigh- Jeans law. 

IRAC data, when combined with the optical, turn out 
to be a powerful tool to identify AGNs, normal galax- 
ies, and stars. In the case of AGNs, the near-IR emis- 
sion is dominated by a dusty torus heated by the cen- 
tral engine, while for normal galaxies, IRAC detects star 
light. Consequently, type-1 AGNs have flatter optical-NIR 
slopes and bluer colors; these kind of sources occupy a 
small characteristic locus in color space, while galaxy col- 
ors vary significantly as a function of redshift. The optical- 
IRAC color-color diagrams are only partially effective in 
disentangling galaxies and stars. In fact, stars and ellipti- 
cals tend to have similar colors at the B — R reddest end, 
i.e. for high redshift ellipticals and cool stars, because the 
spectral continuum in early-type galaxies is dominated by 
old stellar populations. 

Those color-color diagrams that include MIPS 24 /im 
data would be potentially very useful to further distin- 
guish between different object classes, but SWIRE imag- 
ing is not deep enough to detect many normal galaxies at 
24 /xm. This is particularly true for elliptical galaxies. 

6.3. Ultraviolet sources 

Ultraviolet light is mainly emitted by young hot stars 
in galaxies, especially during star formation events. 
Nonetheless, powerful starbursts are typically hosted by 
dusty, thick environments (molecular clouds); in these 
cases, the UV emission of young stars is heavily extin- 
guished and reprocessed to the mid- and far-IR. 

The spectral energy distribution of type-1 AGNs con- 
tains a significant feature in the far-ultraviolet to optical 
region, known as "the big blue bump" . This feature is usu- 
ally explained as the thermal emission from the optically 
thic k accretion disk feeding the central mas sive black hole 
fe.g. IShieldslll978l: iMalkan fc Sargentlll9 82'>. which should 
pea k in the extreme ultraviole t (EUV, 100 < A < 912 A, 
e.g. iMathews k, Ferlandlll987|) . 

The GALEX deep catalog in the ESIS area contains 
~19000 near-UV (1771-2831 A) and ~9000 far-UV (1344- 
1786 A) sources. Roughly 7500 objects are detected both 
in the NUV and FUV channels. GALEX flux densities are 
typically brighter than ~1 [/iJy] in both bands. 

As a result of the cross-correlation between multiwave- 
length catalogs described above, 856 sources turn out to 
be detected both by the NUV GALEX survey and the 
MIPS 24 ^m channel; 803 of these have a secure optical 
counterpart, detected in all the three B, V, and R ESIS 
bands. 

The left top panel in Fig. El represents the distribu- 
tion of sources in the [B — R) vs. {NUV — B) color space. 
Overlaid are the same template tracks as in Fig. 1151 trun- 
cated at a redshift of 1.5. Galaxies and stars are prop- 



erly distinguished, the latter having {NUV — B) > 2 and 
B ~ R < 1.5. The Type-1 AGN template occupy a very 
small area in color space, near to {NUV — B) — {B — R) — 
0. For redshifts smaller than 1, the star burst template 
(dashed line) traces the area with 2 < {NUV — B) < 4 and 
increasing {B — R) > 1; at larger redshift, the {B — R) color 
decreases again and the M82-like templates intersects the 
locus of z < 1 spiral galaxies. 

Apart from stars, the bulk of observed sources lies 
blueward of normal spiral galaxies. These are likely blue 
galaxies detected by GALEX thanks to some moderately 
enhanced, unobscured star formation, or type-1 AGNs. 
The locus of z < 1 starbursts is not very populated, since 
dust extinction allows only the brightest sources to be de- 
tected in the UV. 

The thicker points in the left top panel of Fig. El rep- 
resent sources detected not only in the UV, but also at 
24 /im by Spitzcr. The right top panel of Fig. 1161 com- 
pares the NUV and mid-IR excesses for these objects and 
include also a few stars (at bottom). 

It is worth noticing that, once more, the introduction 
of 24 data is critical in disentangling different kind 
of sources. Template tracks nicely trace the increase of 
ongoing star formation from normal spirals to starbursts. 

The bottom panel of Fig. El shows the ratio of mid- 
IR and far-UV fluxes as a function of the UV color of 
GALEX-24/im sources. The M82 starburst template lies 
above {FUV - 24) > 8: GALEX is not sensitive to young 
stars enshrouded in thick dusty clouds, which instead are 
detected in the mid-IR domain, thanks to dust reprocess- 
ing of the extinguished light. Blue UV galaxies with some 
24 /im excess emission populate the upper part of the dia- 
gram (reddest [FUV — 24]), in between the normal spirals 
and starbursts loci. 

These color-color diagrams suggest that the majority 
NUV-I-MIPS sources are low or moderate redshift late- 
type galaxies {z < 0.5), with enhanced star formation 
activity. This ongoing star formation is almost unextin- 
guished, or moderately absorbed {Ay < 0.5), since both 
an ultraviolet and a mid-infrared excesses are detected in 
the SEDs of these galaxies. 

6.4. SEDs examples 

Multiwavelength coverage from the X-rays to the far-IR 
allows the study of the nature of individual sources. 

Blue optical restframe emission is dominated by young 
stars, while redwards (V, R bands), moderate-aged (type 
A-F) stars are preferentially detected. 

Near-IR imaging traces the old stellar population of 
normal galaxies, and at A > 4 /im, hot (T ~ 1000 [K]) 
dust emission, if present, is no longer negligible (e.g. in a 
AGN- heated torus). 

Mid-IR radiation (5-30 /xm) is dominated by warm 
(T ^ 100 — 300 [K]) dust, either heated by young stars 
in starburst galaxies or belonging to AGN tori. Polycyclic 
Aromatic Hydrocarbons molecules characterize mid-IR 
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Fig. 16. Color-color plots for GALEX ultraviolet sources in the ESIS current area. Template tracks are the same as 
in Fig. El The thick dots in the left top panel represent sources dete cted in the 24 /xm Spitzer channel. The reddening 
arrow refers to the standard extinction law bv lCardein et all lll98flt) . 



restframe light: bright emission bands are detected at 3.3, 
6.2, 7.7, 8.6 and 11.3 /xm and are typical of starburst galax- 
ies. 

In the far-IR spectral domain (30-200 /im), warm 
(T - 100 - 300 [K]) and cold (T - 40 [K], belonging 
to cirrus) dust contribute to the energy budget, produc- 
ing very bright luminosities in active galaxies (starbursts, 
AGNs) as well as in normal spirals. 

Ultraviolet emission comes from young stars, but is 
extinguished by dust and reprocessed to the infrared (see 
also previous Section) . Elliptical galaxies produce a signif- 
icant "UV excess" , with regards to old stellar continuum, 
due to the presence of evolved components such as plan- 
etary nebulae a nd AGB manque stars (see, for example, 
lO'ConnelllHool . 

In the X-rays, AGN power-law emission is detected, ei- 
ther direct (type-I sources) or partially absorbed (type-II) 
in addition to scattered light (if the optical depth is high). 



Star burst gala xies hos t seve ral X-ray components (see, 
e.g., IPersic fc R .cphacli l2002l) . most prominently thermal 
emission by X-ray binaries. Hot plasma in galaxy clusters 
produces extended thermal bremsstrahlung X-ray radia- 
tion. 

Finally, adding deep radio observations will provide 
further insights into the nature of the ESIS galaxy pop- 
ulations, detecting synchrotron radiation from starbursts 
and radio loud AGN. 

Figures El and UHl show the SEDs of five ESIS nearby 
sources, exemplifying different classes of objects. None of 
these sources were detected by ISO, nor by SWIRE's 70 
and 160 /im observations. 

Source no. 1477 is a typical spiral galaxy with some 
enhanced, absorbed, ongoing star formation, displaying a 
moderate mid-IR excess in Spitzer's 8.0 and 24 /im chan- 
nels, and red UV-optical colors. Overlaid is a M51 spiral 
template shifted to z = 0.5. 
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Fig. 17. Spectral energy distributions of three sample sources, from the UV to the mid-IR. The superimposed templates 
are: 1477, spiral galaxy at z = 0.5; 7160, starburst (dot-dashed line) and blue spiral (dashed) at z — 0.8, 76348 spiral 
galaxy a.t z — 0.18. 



Source no. 7160's SED resembles a blue starburst, 
showing a bright IR excess, likely powered by warm dust 
and PAH features, and moderately blue optical colors. The 
overlaid dashed line is a blue spiral template based on the 
nearby galaxy NGC4490 (PoUetta et al., in prep.) shifted 
to z = 0.8, while the dot-dashed line represents a M82 



template at the same redshift, which reproduces the ob- 
served IR excess fairly well. The M82 template emission 
seems to be too absorbed in the optical-UV, being sys- 
tematically fainter than the observed fluxes. This effect 
indicates that this source probably hosts a moderately ex- 
tinguished starburst. 
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Fig. 18. Spectral energy distributions of two X-ray detected sources. Object 81115 has a typical type-1 AGN SED; 
the optical emission from source 81218 is reproduced by an elliptical galaxy template dX z — 0.2 (dashed line). 



Object no. 76348 is a normal spiral galaxy, with some 
weak unextinguished ongoing star formation producing 
some UV excess. This source lies in the area observed with 
SOFI in the J and Ks bands. The dashed line is a MlOO 
template, redshifted to z = 0.18. 



FigurelT^shows two objects detected in the X-rays and 
in the near-IR: a type-I AGN (81115) and a spheroidal 
galaxy (81218). The former SED is a flat power-law spec- 
trum (A-F'(A) = const) over the entire wavelength range 
considered. The latter source shows some 24 /^m emis- 
sion in excess to what is predicted for a simple old stellar 
population (the dashed line represents an elliptical galaxy 
template at z = 0.2), possibly powered by a hidden AGN, 
which is responsible also for the X-ray flux. Note that 
the presumably high UV-optical extinction of the possible 
type-2 AGN is confirmed by the quite hard X-ray spec- 
trum. These plots report \F[X) instead of flux density, in 
order to better highlight the energy budget in the X-rays, 
compared to the other wavelengths. 



7. Summary 

Optical B, V, R imaging of 1.5 square degree in the 
ELAIS-Sl area, belonging to the ESIS survey has been 
presented. The data consist of deep Wide Field Imager 
(2.2m ESO/MPI telescope in La Silla) observations, re- 
quiring >30 different frames per WFI pointing, for a total 
exposure time of 2.5h per band. Ad hoc technical solu- 
tions have been developed for correcting illumination gra- 
dients by means of super- sky fiat fields, and to remove 
color effects and astrometric distortions. Accurate simula- 
tions have been run in order to estimate source extraction 
efficiency and photometric uncertainties. 

A total of 132712 sources are detected on the 1.5 [deg^] 
area analyzed here. Equatorial coordinates of objects are 
constrained within a ^0.15 [arcsec] r.m.s. uncertainty, 
with respect to the GSC 2.2 catalog; flux errors are es- 
timated to be of the order of ~2, 10, 20% at mag. 20, 23, 
24; we reach 95% completeness B,V~25 and R'^24.5. 

Number counts in the B and R bands were compared to 
literature data and theoretical models of galaxy evolution. 
Our data are consistent with previous w orks from vari- 
ous surveys (e.g. HDFN, HDFS, VIRMOS. IWifliams et alJ 
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I199(t1 iMetcalfe et aljEioil: iMcCracken et al.ll2n0.'^ . The 

stellar contribution to the counts has been taken into ac- 
count both using S Extractor's stellarity index and the 
Ijarrett et al.l l)l994j) Milky Way m odel. Comparison to 
number count models by IMetcalfe et al.l |200lll. semi- 
analy tical ACDM mock catalogs f GallCS jHatton et al.l 
1200,11) and ACDM (Nagashima ct al. 200j) predictions 
confirm the need of galaxy luminosity evolution to prop- 
erly reproduce the observed data. ESIS galaxy number 
counts suggest that at z > 1 the rate of galaxy encounters 
was probably larger than today, and the mean restframe 
UV-Blue luminosity of galaxies was larger than locally, 
because of enhanced star formation activity, possibly trig- 
gered by mergers. 

The ESO-Spitzer Imaging extragalactic Survey is pro- 
viding optical identifications, colors, rough morphologies, 
photometric redshifts, for a good fraction of the ^ 300000 
IR sources detected by SWIRE/Spitzer (A =3.6, 4.5, 5.8. 
8.0, 24, 70, 160 Aim) in the 5 [deg^] ELAIS-Sl area. 

In addition to the Spitzer/SWIRE survey, ancillary 
data include observations in the X-r a.vs (XMM, Puccett i 
et al., in prep.), ultra-violet fCALEX. lMartin et alJ2005l) . 
near-IR (NTT, Dias et al.. Buttery e t al., in prep.), and 
radio (ATCA, ICruppioni eraD Il999|) spectral domains. 
Some of the enormous potential stored in this multi- 
wavelength dataset has been traced in the last Section 
of this paper. 

Thanks to the wide wavelength baseline covered, 
Spitzer-optical color-color diagrams turn out to be very 
powerful tools to disentangle different source populations, 
such as normal spirals, elliptical galaxies, starbursts and 
active galactic nuclei. 

The ultraviolet to mid-IR color space was analyzed, 
matching ESIS sources to the GALEX Deep Imaging 
Survey catalog released in January 2005. Roughly 80% 
of the UV sources in the ESIS area are detected in the 
optical, while only a small fraction have a mid-IR 24 /im 
Spitzer counterpart. UV-optical and UV-24/im color-color 
diagrams suggest that sources detected both by GALEX 
in the ultraviolet and by MIPS at 24 /im are low or mod- 
erate redshift late-type galaxies, with enhanced, almost 
unextinguished ongoing star formation activity. 

Finally, X-ray to far-IR spectral energy distributions of 
some ESIS sources were presented as examples of the char- 
acterization of SWIRE objects, based on ancillary data. 

Such a comprehensive wavelength coverage of a large 
area will give a significant contribution, among others, 
in: (1) deriving the physical properties (e.g. ongoing rate 
of star formation and assembled stellar mass) of Spitzer 
galaxies up to 2; ~ 3, (2) searching for distant (z > 1) 
galaxy clusters, (3) studying the statistical properties (e.g. 
number counts) of starburst, evolved galaxies, AGNs and 
comparing them to theoretical models, (4-) studying the 
cosmic star formation density as a function of redshift, 
(5) building the luminosity and mass functions of galaxies 
and studying their dependence on redshift and environ- 
ment. 



Updates on the ESIS project are available on the web 
page http : / / dipastro .pd . astro . it/esis| The data de- 
scribed in this paper are included in the third SWIRE 
Spitzer Legacy data release (Fall 2005). 
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